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1 Introduction 

t> 

The study of substructuring properties of clusters of galaxies is one of the 
basic problems on the way of understanding of the mechanisms of formation 



> 

■^j- ; of the large scale structure of the Universe. It is determined by the fact that, 

the dynamical time scales of clusters of galaxies are comparable with their 
ages, and hence their hierarchical properties do contain essential primordial 
' information. 

A number of statistical methods are developed for the investigation of 
the substructures of galaxy clusters jR], |(| gj. Most of the methods are 
based on the analysis of the positional information, combined with the use 
of redshift information. For example, wavelets (see e.g. [Q] ) are usually used 
to analyze the data on 2D coordinates of galaxies, though the wavelets can 
£2 \ be applied for ID redshift information as well. However such procedure is 

not self-consistent and can lead to various biases. 

From this point of view the S-tree technique developed by Gurzadyan, 
Harutyunyan and Kocharyan (see|5|, |(|) is using in a self-consistent way 
$— i ' the positional, kinematical and magnitude information about the system 

by means of consideration of the phase space of the N-body systems and 
revealing its typical properties. S-tree method has been applied already for 
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the study of substructure of the Local Group 0, the core of Virgo cluster [11], 
and Abell clusters of ESO Key Program survey (ENACS)||. 

The existing version of S-tree || is based on the 'point' approximation of 
the A-body system, i.e. it neglects the dimensions of galaxies. However, as 
distinct to stellar dynamics, when the point approximation is usually quite 
a fair one, the 'non-point' features of galaxies can be ocassionally important 
in the study of clusters of galaxies. 

Below we will inquire into the development of S-tree technique, in non- 
point approximation, namely, considering the galaxies as spheres of given 
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radius, so that the cluster should appear as a set of N spheres. 

Our second aim is the description of a refined algorithm of analysis of 
the numerical output information of S-tree. 

2 S-tree technique 

First, let us briefly summarize the basics of S-tree technique; for details we 
refer to the original papers |], [l]] and to earlier references therein. 

Consider iV-body system as X set. X = {xi, X2, %n}, where Mxi, (i = 
1, . . . , N) are given as follows 

Li(lj,lf,lf) Li is coordinates vector; 

Viiy^i Vi , Vi 3 ) Vi is velocity vector 

and Mi are the masses of particles. 

Consider the following P function : 

P:IxX-t R + and p € -R+. Now we give the basic definitions. 
Definition 1.1 

We will say, that VxGl and My € X are p-bound, if: 

P(x,y) > p. 

Definition 1.2 

We will say, that U C X(U ^ 0) is p-bound group, if: 

a. Vief/ and My G X \U =>- P(x, y) < p 

b. Vx e U and My £ U 3x = Xi x , Xj 2 , . . . , Xi k = y, 
such P(xu, xu +1 ) > p Ml = 1, . . . , k - 1. 
Definition 1.3 

We will say, that U\,...,Ud is the distribution of set X by p-bound 
groups, if: 

a. Uti Ui = X 

b. i^j (i = l,...,d;j = l,...,d)=>U i nU j = Q 

c. MU c i (i = 1, ■ ■ ■ , d) is p-bounded group. 

Note that any P function can be represented as A(aij) matrix Q. 

Thus, we can have a representation of any A^-body system via S'-tree 
diagram algorithm for given A matrix and Mp. This distribution will satisfy 
the definitions 1.1,1.2 and 1.3. This algorithm is efficient for any given p, 
so that one will obtain the evolution of dependence of groups on p. We will 
represent the final results in the form of tree-graph (S-tree) or a table. The 
way one has to choose the final distribution is described below. 
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3 Distribution via D matrix 



Consider again the set X. Via Mj let us denote the projection of Mj on l l l 2 
plane. Introduce the following set: Q = {Q\(M\, Ri), . . . , Qn(Mn, Rn)}, 

where Ri G R+, i = 1, . . . , N and 

Q t (M t ,R t ) = - l]f + (/ 2 - l\f < P 2 }. 

Consider P function, as: 

P.X^Q, P{x i ) = Q i {M i ,R i ) 

For Pp : Q x Q — > i? + and V/? G i2+ we will give definitions, which are 
equivalent to definitions 1.1-1.3. 

Definition 2.1 We will say, that Vxj <G X and Vxj <G X are p-bound, if: 

Pp(Qi,Qj)>P, (i^j) 

Definition 2.2 We will say, that U C X(U ^ 0) is p-bound group, if: 
a. Vxi G U, Mxj G f? Pp(Qj, Qj) < p 

6. \/x,i G f/,Vxj eU(i ^ j) 3xi = Xi ,x h ,. ..,x ik = xj,Vt G (0, . . . fe-1) 

Pp(Qit,Qu+i) > p 

Definition 2.3 We will say, that U\, . . . , is the distribution of the 
p-bound groups of the set X, if: 

a. Uti Ui = X 

b. U i nU j = V),iii^j,Vi,j = l,...,d 

c. VC/j i = 1, . . . , d is p-bounded group. 
Note, that Pp(Qi, Qj) = Vi. 

Obviously, any Pp function can be represented as D(dij) matrix. Now 
we will show a way of construction of the matrix D. 
Ui = j Vi,j = l,...,N dij = 0. 

d . f o if Qi n = 

lj \ s Q 1 nQ J otherwise 

where Qi = Qi\ Tj, Tj-is bound of Qi and Sq^Qj is the area of intersection 
of Qi and Qj. When we know Ri, Rj, Mi, Mj values, it is easy to calculate 
the area Sq^q^ 
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So, using S*-tree diagrams algorithm, we can obtain the distribution by 
p-bounded groups of X set according to defintions 2.1-2.3. 
If we introduce also the following function 

H = Pp o P, 

where P : X x X ^ Q x Q. 

P(xi,Xj) = (p(xi),p(xj)) = (Qi,Qj), then we can use only 1.1-1.3 defi- 
nitions, so long as 

H : X x X -> R 



4 Output information algorithm 

The results of 5-tree analysis can be represented in a form of a table ||, [|. 
For optimization and simplification of the search of final distribution of the 
system, i.e. the search of the necessary floor of the table, one can use the 
following algorithm. 

First, one has to perform a proper transformation of the latter table. 
This will enable us to obtain the information on the quantitative structure of 
each distribution in a convenient form. This transition needs M xN actions, 
where (M x N) is (oy) matrix dimension. We correspond a distribution of 
Ui numbers to each floor, so that 

N 

U i = ^2 a ij> i = l,...,M. 

It is easy to see, that the first and the last floors of the matrix (a^) are 

Ux = N 

N 
3=1 

In some sense Ui is reflecting the distribution 'density' of the system's sub- 
groups. In order to decrease the time of the run of algorithm it is convenient 
to stop the investigation of matrix at some U p , where p is the number of the 
floor with Up > S. One can choose, for example, as x of Um, so that x will 
depend on the matrix dimension. 
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By the next step we construct a chain of auxiliary matrices, in order 
to obtain the final K matrix reflecting the dynamics of quantitative and 
qualitative distribution of changes. 

Each ith matrix is reflecting the quantitative dynamics of ith subgroup. 
At the same time, ith matrix can be represented as a function: 
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where i = 1, . . . , N, t = 1, . . . , M and n\ is the number of zth subgroup at 
ith floor. 

Note, that if we 'paste together' all mentioned auxiliary matrices by 
'floors', we will obtain exactly the sought-for K matrix. 

Now let us turn to the matrix K. We will consider the behavior of the 
subgroup, which contains maximal points (m > 2). We will fix Api interval, 
i.e. A/jj-number of floor, when the group is preserving its main structure or 
if the group is loosing, for example, the x structure. Note, that x has to be 
defined for each matrix: if the group is loosing more than x of its numbers, 
then the process will continue up to the next floor. 

Obviously, this process is finite, so far as the number of subgroups Cj is 
limited: 

ci < a < cm 

, where c\ = 1 and cm = N, i.e. the first floor contains all particles from a 
single group (whole system), while on the last floor each particle represents 
a separate group. This readily follows from the Definitions 1.1-1.3 and the 
method of choice of pt- Then, we have to choose the maximal Api and will 
take the corresponding distribution. As a final distribution we will consider 
the floor, which results after the destruction of the mentioned distribution. 

For the realization of this algorithm we need about D = M a x N b actions 
||, where a, b > 0. 

Let us schematically illustrate this algorithm on a concrete example. 
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5 Conclusion. 

We have proposed some new ways of development of S-tree method in order 
to take into account the effects which could be important for description 
of clusters of galaxies. Namely, we considered the possibility of accounting 
for the non-zero dimensions of the galaxies. The described generalization 
of S-tree can enable the construction of new type N-body models and can 
be important especially while studying the denser core regions of galaxy 
clusters. 

We also represented an efficient procedure for the reliable analysis of 
the output information of S-tree. As numerical experiments, as well as the 
S-tree runs of the real galaxy clusters indicate, this procedure can become 
especially useful, when the numbers of galaxies in the studied samples can 
become about 1000 and more. The obtaining of samples over thousand 
of galaxies with known redshifts and magnitudes in the vicinity of a given 
filament, is within the aims of several forthcoming observational programs 
(for reviews see Durret et al 1994). 

We thank V.G.Gurzadyan and K.Oganessyan for valuable discussions. 
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